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AB5TRACT 



An investigation of using the response of a circuit at. selected test 

frequencies to isolate faulty circuit components it made.. A procedure 

using a sensitivity approach for intelligent selection of test frequencies 

is developed. The developed procedure is tested anct the results are 

compared with results using conventional procedures!.. The effect of 

» 

random, within tolerance variations of nonfaulty components on the 
results is studied for both conventional and developed procedures . 
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I. INTRODUCTION 



Maintenance of electronic equipment is a problem that increases as 
equipment becomes more complex. With ever increasing costs (of equip- 
ment) it also becomes imperative to decrease equipment downtime as 
much as possible. To this end a fast, reliable method of locating 
circuit failures must be used. 

"It has been customary to isolate failures by testing individual compo- 
nents. This is usually done by isolating the fault to a particular stage, 
followed by subassembly isolation and finally component isolation by 
individual testing. 

With the availability and flexibility of digital computers, a pos- 
sible solution to the problem of excessive downtime is to use a computer 
to help isolate the fault. The first step in accomplishing this is the 
development of a testing and isolation procedure which could be program- 
med for the digital computer. 

A recently presented procedure has been examined by a number of 
individuals [1], [2], [3], and [4]. The procedure was described by Seshu 
and Waxman [1] with other investigators testing and making applications of 
the procedure with few modifications. 

The basic procedure is to examine the response of a circuit under test 
at selected frequencies. Comparisons of the results with a list of "library 
results" can pinpoint a faulty component. 

Most investigators have used test frequencies near the poles and zeros 
of the function under consideration. Other frequencies have been chosen 
between the pole and zero frequencies [2] and [4]. In this study an 
alternate method of test frequency selection was developed and tested and 
results were compared with results using other methods. 
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The procedures for signature comparison used by other investigators 
have not been explained in detail in most reports. However, it appears 
that a nearest neighbor approach is the most common technique used to 
date. Two alternate methods were examined and are discussed with 
conclusions as to their practicality. 

Previous methods of fault determination after signature comparison 
have also not been documented. Two possible methods were tried. The 
first method used computer selection for each network function and manual 
correlation for fault isolation. The second method was completely program- 
med for digital computation. 

As far as can be determined from available reports no investigator 
approached the problem of error due to in-tolerance variations of com- 
ponents. This study included an investigation of that problem with both 
the conventional and the proposed procedures. 



10 



II. PROCEDURAL BACKGROUND 



A. FUNDAMENTAL CONCEPT 

All components of a given network play a part in the formation of any 
of the network functions. The sensitivity of any given network function 
to a particular component is generally unique to that component. That is, 
in a network with K components a network function will have K sensitivity 
functions, one for each component. A change in the value of a particular 
component may cause a decrease in the magnitude of a network function at 
one frequency and an increase in the magnitude at a different frequency. 
There may in fact be frequencies at which the function sensitivity is 
near or equal to zero for a particular component, excluding catastrophic 
failures. 

Specific knowledge about this behavior can be used to detect and 
isolate component value changes. If, for example, it is known that the 
input impedance is insensitive at low frequencies to changes in a specific 
component, one would assume that component was not faulty if the input 
impedance exhibits a significant change at low test frequencies. 

B. POLE-SHIFT TECHNIQUE 

The component variation effect on frequency response can be observed 

with the pole-shift technique. A network function can be represented by 

a constant and a unique set of poles and zeros; i.e., 

fM . y ( S -V(S-Z 2 ) (S-V 
F(S) ' K (s- Pl )U-p 2 ) ... (s-p n ) 

where 

F(s) = any network function 
K = constant 
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z.j = zero frequency 

Pi = pole frequency. 

Every component plays a part in pole and zero locations. However,, a 
few components may dominate the location of a given pole or zero frequency. 
Thus, if one could detect the shift of the poles and zeros of a gfven 
function, it may be possible to determine component variations. A simple 
example can best illustrate this point. 

Example 1 - Consider the network shown in Figure 1. 




Figure 1. Example 1 - Low Pass Filter 
Vn 

The network function -r— (s) is given by: 

A 1 



Zi= s 2 + 11s + 11 
*1 s 2 + 10s + 1 

_ (s + 9.887) (s + 1.113) 

(s + 9.899) (s + 0.101) * 

The pole and zero shifts were determined by varying each component 
value by ten percent and observing the corresponding critical frequency 
shifts. The results: 
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Circuit 

Condition 


h 


Location 

P 2 


(radians/sec) 
Z 1 Z 2 


Pi 


Percent Shift 
P 2 Z T 


Z 2 


Nomi nal 


.101 


9.899 


1.113 


9.887 


0 


0 


0 


0 


R-j variation 


.101 


9.899 


1.020 


9.889 


0 


0 


.61 


.02 


C-j variation 


.092 


9.908 


1.010 


9.899 


8.9 


.09 


9.15 


.12 


L.j variation 


.101 


8.990 


1.114 


8.977 


0 


9.2 


.08 


9.2 


R2 variation 


.092 


10.908 


1.101 


10.899 


8.9 


10.2 


T.07 


TO. 2 



TABLE 1 1 -1 . Percent Pole and Zero Shifts 



It is immediately obvious that a change in or C-j will cause a greater 
shift in Z-| than will a corresponding change in L-j or . Also, p2 and Z2 
are most affected by and R2 and p^ is affected only by C-j and R2. 



C. CHANGES IN RESPONSE FROM POLE SHIFT 

In practice the exact location of the critical frequencies can usually 
not be determined directly from the circuit. However, this does not rule 
out the pole shift technique. Any shift in critical frequency causes a 
change in the magnitude of the given function at all frequencies greater 
than the shifted one. This can be seen from the simple Bode representation 
of a pole: 
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If f is changed to f' (indicated by dashed line), the function 
c c 

magnitude, at any frequency, f , greater than f , is changed (increased 
in this case). The magnitude is not affected at frequencies less than 



Now, from the network in Figure 1 we obtain: 



V 

2 _ = 10 = 10 

^1 s 2 + 11s + 11 (s + 9.887) (s + 1.113) 



Figure 3a is the Bode plot for the function i^/Vi . Figure 3b and 
3c show the effect of increasing and of increasing , respectively . 
The nominal response is shown as a dashed line. 
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It can be seen that an increase in R-j causes a decrease in the 
magnitude of the response between and p^. A lesser decrease is seen 
at frequencies greater than P2« The response between p-j and P2 is 
increased for an increase in L-j while the response above P2 is decreased. 
Thus, if one examines the response at a frequency between p^ and P2 and 
at a frequency greater than Pg it will be a reasonably simple task to 
distinguish between an increase in R-| or L-j . This example is an over- 
simplification of the procedure, but the application of which can be 
appreciated. 
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III. GENERAL PROCEDURE 



A. BASIC PROCEDURE 

The basic procedures outlined by Seshu and Waxman [1] was followed 
with modifications as explained herein. Two network functions were used 
concurrently. Work done by Maenpaa, Stehman, and Stahl [5] indicated 
that the redundancy from using two functions was more than adequately 
compensated for by the increased fault detection capability. 

The symbolic network functions were calculated by hand. It should 
be pointed out that computer programs exist which calculate symbolic net- 
work functions. These programs should be used for more complex circuits. 

The functions were numerically evaluated and factored to determine 
the poles and zeros. Test frequencies were then chosen. Previous 
procedures generally selected test frequencies based on the pole and 
zero locations. As a test of the procedure several test frequencies were 
chosen between the critical frequencies as well as 0 hertz and one 
frequency above the highest critical frequency. The actual selection 
was programmed for computer selection and is described in Section VII. 

It was felt that the procedure of selecting test frequencies based 
on poles and zeros was, although intuitively satisfying, a somewhat 
illogical procedure. The concern was not with the function itself but 
with the behavior of that function with a parametric change. Thus, 
sensitivity functions were calculated for each element. Based on these 
functions, six test frequencies were chosen. This procedure is discussed 
in Section IV. 

After selection of test frequencies, a worst case response, for 
component values within 10 percent tolerance limits, was computed for 
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each test frequency. This defined a nominal range of performance. This 
is not necessarily the best method for definition of nominal range. This 
range might be specified for a particular circuit. However, in the 
absence of such specifications the above procedure was followed. 

Each component was then allowed to take on several discrete values 
out of tolerance limits. All other components were held fixed at their 
nominal value. The response at each test frequency was then calculated 
and catalogued for the entire collection of sets of component values. 

The catalogued response values were then quantized into nine levels 
and labeled one through nine. Quantization level five was reserved for 
response within nominal limits. The setting' of the levels is explained 
in Part B. The quantized responses for a given set of component values 
is defined as the signature for the given circuit configuration. 

The signatures were computed for the following cases: 

value = zero 
+ 20 % 

+ 50 % 

100 times nominal. 

The entire set of signatures for a given function is defined as the 
library for that function. There is no theoretical limit, on the size 
of the library. Hov/ever, the inclusion of many more circuit configurations 
would probably only complicate the matching procedure and not increase 
detection capability. 

Using the primary test circuit the input impedance and voltage gain 
at each test frequency was measured. Then, after quantization, the signa- 
tures measured were compared with the library signatures and a "closest 
match" was made. Thus, the faulty component was identified. 

In order to test the procedure, several error conditions were set and 
the response calculated. In order to more closely approximate actual 
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circuit conditions, all component values were allowed to vary randomly 
(uniform distribution) within 10 percent tolerance limits. Calculations 
were made with and without the random variations to determine how much 
effect this would have on the results. 

The results of all tests are presented in Section VI. 

B. SETTING OF QUANTIZATION LEVELS 

. With the exception of selection of test frequencies the most critical 
step in the procedure was the setting of quantization levels. Two methods 
were used successfully. One was programmed as part of the main analysis 
program. The second method, although not programmed is possibly adaptable 
to computer selection. Both methods have nine levels of response with 
level five reserved for the nominal range. 

•The nominal range was selected by taking the range between the 
maximum positive and maximum negative deviation from normal for any 
element value at j^lO% from normal. The programmed method then set the 
other eight levels as shown in Table III-l in which F is a multiplying 
factor used to force the possible response range to cover all quantum 
levels with levels one and nine assigned only for extreme cases. 

Normalized Response 

(-,1+13 FN] 

(1 + 13 FN, 1 + 7 FN] 

(1 + 7 FN, 1 + 4 FN] 

(1 + 4 FN, 1 + N] 

(1 + N, 1 - M) 

[1 - M, 1 - 4-FM) 

[1-4 FM, 1-7 FM) 

[1 - 7 FM, 1 - 13 FM) 

[1 - 13 FM, 0) 

Quantization levels. 
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Quantum Level 

1 

2 

3 

4 

5 

6 

7 

8 
9 

TABLE III-l. 



The second method used a more inituitive approach. The maximum 
deviation levels, 1 and 9, were selected based on the maximum deviation 
from normal under any condition. The mid levels were then selected so 
that the catalogued or library signatures would cover all quantum levels 
as uniformly as possible. 

The manual approach is acceptable for small circuits but gets 
unmanageable quickly as circuit size increases. Parts of the procedure 
could be easily programmed. However, a suitable algorithm would have to 
be developed prior to complete programming. 

C. TWO COMPONENT VARIATIONS 

It cannot be assumed that only one component at a time will undergo 
changes in value, "in order to investigate the possibility of tv/o compo- 
nent values changing simultaneously, two possibilities were considered. 

If one component masks function changes due to variations in another 
component (i.e., no change in function magnitude), the entire procedure 
will not work. However, if one is primarily concerned with circuit 
performance, the failure to detect component variations is of no concern 
if there is no degradation in circuit performance. 

The second situation considered was one in which there was an obvious 
change in circuit response due to two components varying. The original 
procedure made allowances for varying only one component at a time. How- 
ever, a slight modification allowed for the generation of a set of library 
signatures for two component variations. The set is presented in the 
computer output section. 

It soon became apparent that in order to use the two component 
variation library, one needed apriori knowledge of the number of faulty 
components. Since this is not usually available, the usefulness of two 
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component signatures is doubtful. Additionally, an examination of the 
library signatures revealed that one of the two components varied usually 
dominates the response. Thus, if two components have failed the single 
component library would probably isolate the dominant component. After 
replacement of that component further testing could reveal the second 
faulty component. Consequently no further investigation of multi- 
component variation was undertaken. 
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IV. PRIMARY TEST CIRCUIT 



The primary circuit used throughout the investigation was a srx- 
element low-pass filter. (See Figure 4.) The unsealed frequency 
response, represented by the voltage ratio transfer function,. V^/Vy, 
is shown in Figure 5 and the unsealed input impedance as a function - 
of frequency is shown in Figure 6. On Figures 5 and 6 the poles are 
-represented by X and the zeros by 0. 




The component values were magnitude scaled by 100 and frequency 

5 

scaled by 10 . This yielded component values of approximately the 
same magnitude. Ideally, the scaled cutoff frequency shouTd be V Hz.. 
However, with this scaling f = .02 Hz and no difficulty was encounter- 
ed. 

The scaled network functions of interest are: 

V 2 = .1634 

^7 s 4 +l 7 . 67s 3 +6 . 462s 2 +2 . 601 s+. 21 24 

.1634 

(s+17.31 ) (s+. 0997) (s+1 .325+j3.25 ) (s+1 .325-j3.25) 
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w) 



X 



db 



0 - 



- 20 ' 



-40 



.001 



x 



X 




Figure 5. V^/V-j , Primary Test Circuit 
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= s 4 -H7.67s 3 +6.462s 2 +2.601s+.2124 
*1 5.882s 3 +.147s 2 +.817s+.01634 

(s+1 7.31) (s+.0997) (s+1 .325+j3 . 25) (s+1 .325-j3 .25) 
(s+.02) (s+.0025+j .373) (s+.0025-j .373) 

The poles and zeros in Hz are listed in Table IV-1 for convenience. 



V 1 ! 




Vh 


Zeros* 


(Hz) 


Poles (Hz) 


.0518 


(Second Order) 


.0594 (Second Order) 


.0159 




.00319 



2.76 

* Zeros of / 1^ are the poles of V^/V^ . 

TABLE IV-1. Pole and Zero Magnitudes. 
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V. SELECTION OF TEST FREQUENCIES 



A. CONVENTIONAL METHOD 

The conventional procedure was followed with regard to selection of 
test frequencies in order to test the method. It soon became apparent 
that the number of test frequencies was too large to handle with ease. 
Thus, only frequencies between the pole and zero frequencies were 
selected as well as 0 Hertz and one above the maximum critical frequency. 

The computer program written for implementation of the procedure (see 
Program 1) did the frequency selection automatically. The test fre- 
quencies were selected by taking . ,, . 

. f 1 1 - f 1 

f = f 1 + _£ 

'test T c - 2 

• • * i 

for every f^ where f^ = i tn critical frequency. If 



’ . f 1+1 - f 1 ‘ 

f' + c ■ ■ c - 




f i +2 f i +1 

*i+l c " T c 


f c 2 




f c " 2 



one of the frequencies was eliminated. Test frequencies were computed 
for the input impedance because the critical frequencies for V^/Vi are 
included in the critical frequencies for / I ^ . One test frequency was 
eliminated due to the above inequality. This resulted in the selection 
of the following seven test frequencies for a fourth order filter: 



0.0 


Radians 


.001 


Radians 


.00498 


Radians 


.196 


Radians 


.5488 


Radians 


8.49 


Radians 


25.8 


Radians 



TABLE V-l . Conventional Method Test Frequencies 
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B. SENSITIVITY APPROACH TO FREQUENCY SELECTION 

Inspection of the signatures generated using the conventional method 
showed that redundancy existed to the extent that the method was some- 
what inefficient. Two test frequencies, D.C. and the lowest nonzero test 
frequency, gave nearly identical results. This indicated one of the two 
frequencies v/as probably not required. It was also felt that unless 
further investigation was undertaken it would be possible to ignore a 
range of frequencies which could be extremely helpful. 

★ 

This investigation took the form of calculating and plotting functions 
which defined the sensitivity of ^ 2^1 an< ^ *1^1 var i a tions in each 
component. These plots are presented as Figures 7 through 18. Figure 19 
is a superposition of the plots of the sensitivity function for 1 -j/V.j and 
Figure 20 is the same for V 2 /V 1 . Ideally / 1^ should have been used. 
However, the computer program used in this portion of the study did not 
allow for using V-j/I^. It was felt that no loss of information would 
result. 

The superposition plots Figures 19 and 20 are very revealing. The 
first observation is that V 2 /V 1 is not sensitive to changes in R-j and R 2 
at high frequencies nor to C's or L's at low frequencies. . While this 
information is not new the sensitivity functions assist in determining 
when C's and L's have little effect. This can be used to eliminate the 
redundancy observed using the conventional method of frequency selection. 

The mid-frequency range is the most interesting range for the selection 
of test frequencies. For example, near .02 Hz ^ 2 ^] 1S ver ^ sens itive to 
changes in R-j and C 2 , less sensitive to changes in L 2 , even less sensitive 



* The program NASAP-70 was used to calculate the various sensitivity 
functions . 
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to changes in ft, and and virtually insensitive to changes in L-j . 

There is a peaking of the sensitivity functions for , C 2 , and near 
.06 Hz. At this same frequency the sensitivity functions for R-| , R£, 
and L.j are near zero. 

The 1/^ and V-j plots show that at near .3 Hz only the sensitivity 
function for R-j is high. The functions for and are low and the 
other functions are zero. At 3 Hz I-j/V-] is about equally sensitive to 
changes in L. and R-j while the other functions are very close to zero. 

It can also be seen that the frequencies discussed in connection with 
V £ /Vi are reasonably good choices for I ^ /V^ . 

Based on the above observations the following six frequencies were 
selected as test frequencies: 

Number Frequency 

1 .00048 Hz 

2 .027 Hz 

3 .06 Hz 

4 .299 Hz 

5 2.99 Hz 

6 49.97 Hz 

TABLE V-2. Sensitivity Function Test Frequencies. 

After some preliminary work it was observed that the response at .06 
Hz was not well behaved. Since four sensitivity functions peaked for 
I-j/V-j, random variation within tolerance limits of these components were 
being observed as significant changes in response. In order to overcome 
this problem, test frequency three was changed to .055 Hz. The new fre- 
quency was off the peak enough to be useable and still yield the desired 
information. 
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Figure 7. Sensitivity of V^/V-j ^ • 
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• Figure 8. Sensitivity of V 2 /V^ to LI. 
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4. 



|F s (w)| 




.001 .01 .1 1. 

Figure 10. Sensitivity of V 2 /V-| to R2. 



10. w 
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Figure 11. Sensitivity of V ^ /V^ to L2 
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|F s (w) | 




.001 .01 .1 1. 



Figure 12. Sensitivity of V^/V^ to C2. 
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4. 



|F s (w)| 




Figure 13. Sensitivity of 1^ /Vi to R1 . 
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! F s (w) | 
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Figure 15. Sensitivity of I.j/V.j to Cl. 



35 



|F $ (w)| 
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Figure 17. Sensitivity of I-j/V-j to L2. 
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4. |F $ (w)| 



3. 



2 . 



1 . 
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.001 




Figure 18. Sensitivity of to C2. 
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4. | F s (w) | 




Figure 19. I^/V-j Composite Sensitivity Functions. 
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Figure 20. l^/V-j Composite Sensitivity Functions. 



40 



C. COMPARISON OF METHODS 



The comparison of the results using the two methods of frequency 
selection is included in Section VI. Some comments as to appTi'cabil ity 
should be made at this point. 

If one computes sensitivity functions by using the derivative af 
the network function with respect to the element under consideration it 
becomes 



f s u> = 





q(s) 



^ - PM 

k 

[qU )] 2 



dq(s) 

de k 



where: 



f s (s ) = sensitivity function 

d/de. = derivative with respect to element k 



eM 

q(s) 



= ratio of polynomials representing the network function. 



Thus, the poles of the sensitivity function are the same as tha poles of 
the original function. Consequently, if one selects test frequencies on 
the basis of function poles the selection is linked to the sensitivity 
functions. Generally, however, the zeros of the sensitivity function are 
unique and are thus ignored by any selection procedure based on' the poles 
and zeros of the original function. 

From the preceding paragraph, some of the advantages of using sensi- 
tivity functions for test frequency selection become apparent.. However, 
that method has the disadvantage of requiring the calculation erf sensi- 
tivity functions for each element. 
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With the small circuit used as an example the value of sensitivity 
functions are demonstrated. With more complex circuits it appears more 
probable that important test frequencies might be overlooked if the sen- 
sitivity functions are not used. 
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VI. RESULTS WITH PRIMARY CIRCUIT 



A. RESULTS USING BASIC PROCEDURE 

Using the procedures outlined in Section IIT-A and V-A a set of 
signatures was generated for the six selected circuit conditions. 

These signatures are presented in Figure 25. A second set of signatures 
was then generated for the same circuit conditions with random variation 
allowed for the circuit elements. The second set of signatures was input 
to the selection program to determine what errors resulted from the random 
variations. 

Table VI-1 lists the input condition and the program output. The 
numerical result was determined as follows: 

1 = Selection was correct condition 

2 = Selection was correct element 

3 = One of multiple selections was the correct condition 

4 = One of multiple selections was correct element 

5 = Incorrect selection. 

By correlating the results from V 2 /V -1 and V-j/I^ a "hit", "miss", or 
"possible hit" rating was assigned for each test signature as follows. 

A hit was assigned if both V 2 /V 1 and /I^ yielded a 1 or 2 or if one 
was 1 or 2 and the other was 3 or 4. A 3-5, 4-5, or 5-5 combination was 
called a miss. All other combinations were. called possible hits. For 
some cases deviations were warranted. For example, a 3-4 selection could 
have been called a miss if a component other than the given component 
was clearly indicated as being faulty. 

The result was 21 hits, 9 possible hits and 6 misses for 58% hits and 
83% hits or possible hits. Since only the circuit conditions used for 
signature generation were used to generate the test signatures 100% 
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Circuit 


Program 


Numerical 


Program 


Numerical 


Condition 


Condition 


Result 


Condition 


Result 


R1 : Short 


R1 : Short 


1 


Rl : Short 


1 


-50% 


Rl: -50% 


1 


L2: Open 


5 


-20% 


R1 : Open 


2 


L2 : +50% 


5 


+20% 


L2: +20% 


5 


Rl: +20% 


1 


+50% 


Rl: +20% 


2 


Rl: +20% 


2 


Open 


Rl : Open 


1 


L2: Short 
Rl : Open 


3 


LI : Short 


LI : Short 


1 


LI: -50% 


2 


-50% 


LI: -50% 


1 


LI: -50% 


1 


-20% 


LI: -20% 
R2: +20% 


3 


LI: -20% 


1 


+20% 


LI: +20% 


1 


LI: -20% 


2 


+50% 


LI: +20% 


2 


LI: -20% 


2 


Open 


LI : Open 


1 


R2: +20% 


5 


Cl : Open 


Cl : Open 


1 


LI: +50% 


5 


-50% 


Cl: +20% 
C2: -50% 


4 


Cl: -50% 


1 


-20% 


C2: +20% 


5 


Rl: -20% 


5 


+20% 


LI: +20% 
C2 : +20% 


5 


LI: -50% 


5 


+50% 


C2: +20% 


5 


Cl: +20% 
C2: +20% 


3 


Short 


Cl : Short 


1 


Cl : Short 


1 


R2: Short 


R2 : -50% 


2 


R2: Short 


1 


-50% 


R2: -50% 


1 


R2: -50% 
R2 : -20% 


2 


-20% 


R2: -20% 
L2 : +20% 
C2: +20% 


3 


Rl: -20% 


5 


+20% 


C2: +20% 


5 


Cl: +20% 
C2: +20% 


5 


+50% 


R2 : +20% 


2 


R2 : +20% 


2 


Open 


R2: Open 


1 


Cl : Open 


5 



TABLE VI-1. Basic Procedure Results. 



44 



L2: Short 


L2 : 


Short 


1 


L2: 


Short 










Rl: 


Open 


-50% 


L2: 


-50% 


1 


L2: 


-50% 


-20% 


L2: 


-50% 


2 


L2: 


-20% 




LI: 


+20% 




L2: 


+20% 


+20% 


L2: 


+20% 


3 








C2 : 


+20% 








+50% 


L2: 


+20% 


2 


L2: 


Open 


Open 


L2: 


-20% 


2 


L2: 


Open 



C2: Open 


Cl: 


+20% 




R2: 


+20% 


C2: 


-50% 


4 






-50% 


Cl: 


-20% 


5 


Cl: 


+50% 










Cl : 


-20% 










C2: 


-20% 


-20% 


Nominal 


5 


L2: 


+20% 


+20% 


C2: 


+20% 


1 


C2: 


+20% 


+50% 


Nomi nal 


5 


C2: 


+20% 










Cl: 


+20% 


Short 


C2: 


Short 


3 


C2: 


Short 



— R2 : Short 



3 

1 

1 

1 



2 

1 



5 



3 



TABLE VI -1 (Continued) 
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hits or possible hits could have been accomplished. Thus, using the 
conventional procedure, within tolerance variations introduced 17% to 
42% errors. 

B. RESULTS USING SENSITIVITY FUNCTIONS 

Using the frequencies determined from the sensitivity functions a 
set of signatures was generated for the six circuit conditions. In order 
to check the effect of variations within tolerance for this method the 
same procedure as used in VI-A was followed. The result was 50 hits, 16 
possible hits and 6 misses. The 6 misses included 3 due to the test 
signature showing normal response from the circuit. These results showed 
an error rate of 11% - 31%. 

Eight additional circuit conditions were simulated for each element 
to see what could be expected from a less artificial situation. Table 
VI-2. lists all the conditions set and the numerical result. The numerical 
result was assigned as follows: 

1 - closest library signature selected 

2 - right element selected 

3 - one of multiple signature closest to given condition 

4 - one of multiple signatures right element 

5 - incorrect selection. 

Hits, misses and possible hits were also assigned as before. The 
result was 115 hits, 43 possible hits and 10 misses. Or, 68% hits and 
94% hits or possible hits. 

C. USING SUBROUTINE DET FOR CORRELATION 

The same set of test signatures used in VI-B was used to test the 

programmed correlation and selection. The input and results are listed 
in Table VI-3. The results are coded as follows: 
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Circuit 

Condition 


Run 

V 2 /V l 


1 

v^i. 


Run 2 

1 V V 1 


V 1 ! 


R1 : Short 


1 


1 


1 


1 


-90% 


1 


1 


1 


1 


-70% 


1 


1 


1 


1 


-50% 


1 


1 


1 


1 


-40% 


1 


1 


1 


1 


-20% 


1 


5 


3 


1 


+20% 


1 


3 


1 


3 


+30% 


1 


1 


1 


1 


+50% 


1 


2 


1 


1 


+60% 


1 


1 


1 


1 


+90% 


1 


1 


1 


1 


x3 


1 


1 


1 


1 


X75 


1 


1 


1 


1 


Open 


1 


1 


1 


1 


LI : Short 


1 


1 


1 


1 


-90% 


1 


1 


1 


1 


-70% 


1 


1 


2 


2 


-50% 


1 


1 


1 


1 


-40% 


1 


1 


1 


1 


-20% 


1 


1 


1 


1 


+20% 


1 


5 


1 


5 


+30% 


1 


3 


1 


1 


+50% 


1 


1 


1 


1 


+60% 


1 


1 


1 


1 


+90% 


1 


1 


1 


1 


X3 


1 


1 


1 


1 


X75 


1 


1 


1 


1 


Open 


1 


1 


1 


1 


R2: Short 


3 


3 


3 


3 


-90% 


3 


5 


3 


2 


-70% 


3 


5 


3 


5 


-50% 


1 


1 


1 


1 


-40% 


1 


1 


1 


1 


-20% 


1 


1 


1 


1 


+20% 


1 


1 


1 


1 


+30% 


1 


1 


1 


1 


+50% 


1 


1 


1 


1 


+60% 


1 


1 


1 


1 


+90% 


1 


1 


1 


1 


Open 


1 


1 


1 


1 




TABLE VI -2. 


Results 


Using Sensitivity Functions 
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Cl : Open 


1 


1 


1 


-90% 


5 


1 


5 


-70% 


5 


1 


5 


-50% 


3 


1 


.5 


-40% 


3 


1 


3 


-20% 


1 


3 


5 


•+20% 


4 


5 


5 


+30% 


3 


1 


5 


+50% 


5 


1 


5 


+60% 


3 


1 


5 


+90% 


4 


1 


3 


X3 


5 


5 


5 


X75 


1 


1 


1 


Short 


1 


1 


1 



L2: Short 


2 


1 


2 


-90% 


1 


5 


1 


-70% 


1 


1 


1 


-50% 


1 


4 


2 


-40% 


1 


3 


1 


-20% 


2 


5 


3 


+ 20 % 


3 


3 


4 


+30% 


5 


3 


5 


+50% 


5 


1 


3 


' +60% 


5 


1 


5 


+90% 


5 


1 


5 


X3 


5 


1 


5 


X75 


5 


1 


5 


• • Open 


1 


1 


5 



C2: Open 


1 


1 


1 


-90% 


1 


1 


5 


• -70% 


1 


1 


5 


-50% 


3 


1 


4 


-40% 


3 


1 


3 


-20% 


5 


3 


5 


+20% 


5 


3 


5 


+30% 


3 


1 


5 


.+50% 


5 


1 


5 


+60% 


3 


1 


3 


+90% 


3 


1 


5 


X3 


5 


1 


5 


X75 


3 


3 


3 


Short 


3 


3 

TABLE VI -2 


5 

(Continued) 



1 

1 

1 

1 

1 

3 

1 

1 

1 

1 

1 

4 

1 

1 

1 

5 

5 

3 

3 

3 

3 

3 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

3 

3 

3 

1 

1 

1 

1 

1 

3 

3 
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1 - correct primary selection 

2 - one of two primary selections count 

3 - one of two secondary (no primary selections) 

4 - incorrect primary selection. 

The results were 61 rated 1, 12 rated 2, 21 rated 3 and 2 rated 4. 
Or, 63% completely correct and 98% with one of two selections correct.' 
This is roughly the same result using manual correlation. The algorithm 
used for determining primary and secondary selection is outlined in 
Section VII-C. 
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Circuit Condition 



V V 1 V 1 ! 

R1 

-90% 1 1 

-70% 1 1 

-60% 1 1 

+30% 1 1 

+60% 1 1 

+90% 1 1 

X3 1 1 

X75 1 1 



LI 

-90% 1 1 

-70% 1 1 

-60% 1 1 

+30% 2 2 

+60% 1 1 

+90% 1 1 

X3 1 1 

X75 1 1 

Cl 

-90% 3 3 

-70% 1 1 

-60% 3 3 

+30% 3 2 

+60% 3 2 

+90% 1 1 

X3 44 

X75 1 1 



R2 

-90% 1 1 

-70% 1 1 

-60% 1 1 

+30% 1 1 

+60% 1 1 

+90% 1 1 

X3 1 1 

X75 1 1 



TABLE VI-3. Computer Correlation Results. 
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L2 

-90% 

-70% 

-60% 

+30% 

+60% 

+90% 

X3 

X75 



C2 

-90% 

-70% 

-60% 

+30% 

+60% 

+90% 

X3 

X75 

TABLE VI 



3 

1 

1 

2 

2 

3 

3 

3 



3 

2 

3 

1 

1 

1 

3 

2 

3 (continued) 



3 

3 

1 

2 

2 

3 

3 

3 

3 

2 

3 

1 

1 

1 

3 

2 
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VII. COMPUTER PROGRAMS 



Several computer programs were written in Fortran IV to accomplish 
much of the mathematical work associated with this study. Four of the 
programs are included herein with detailed descriptions. Each of these 
programs were used with several modifications. However, only the basic 
program is presented in each instance. All listings are in the designated 
section. 

A. PRIMARY ANALYSIS PROGRAM 

The primary analysis program contains the main program and three 
subroutines. A fourth subroutine, PRQD, is used for root finding and 
is not included in the listing. 

The program calculates the magnitude of both network functions as a 
function of the complex frequency, s, finds the roots of the polynomials 
involved and selects test frequencies. 

Each element in turn is then set to each out-of-tolerance value as 
selected previously, the magnitude of each function at each test fre- 
quency is calculated and all responses are quantized into signatures. 

The output consists of both network functions, the test frequencies, 
the nominal response at each test frequency and the list of library 
signatures. 

The listing is the program which calculates the test frequencies as 
explained in Section IV-A. To modify the program for the selection from 
sensitivity functions the portion of the program between indicated com- 
ment cards was replaced with a read statement and the test frequencies 
were used as the input. 

Figure 21 is the flow graph for the main analysis program. 
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Figure 21. Flowgraph of Primary Analysis Program 
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B. TWO COMPONENT VARIATION PROGRAM 

The main analysis program was modified to generate signatures for 
simultaneous two-component failure. Two different modifications were 
made, one generated signatures for only short and open element failures. 
The other allowed for short, -50%, +50% and open failures. 

Since the programs were only modifications to the main analysis 
program flow data are not included and the listings are omitted. 

C. FAULTY COMPONENT SELECTION 

This program was written to determine possible and probable component 
failures from the two input signatures. The signature from / 1 -j and 
V 2 /V -1 are both processed. A component is flagged as possibly failed if 
one of the library signatures for that component is selected as a match. 

If two or more library signatures for a component are selected as a match 
that component is flagged as probably failed. 

The flow graph for the program is shown in Figures 22-24. The main 
program allows for several signatures to be read in and sent to subroutine 
DET, one pair at a time. 

Subroutine DET first checks to see if the signature indicates the 
circuit is functioning properly (i.e., signature is all 5s). If so, a 
message is printed and control is returned to the main program. The first 
entry into DET with a non-nominal signature results in the library signa- 
tures being read in and stored. 

When it has been determined that non-nominal signatures are being 
processed, both input signatures are compared with the library signatures 
to determine the closest match distance. The signatures are compared a 
second time and those library signatures that are the minimum distance 
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Figure 22. Flowgraph of Main Program for DET 



55 



1 





Figure 23. Flowgraph of Subroutine DET 
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2 




Figure 24. Flowgraph of Subroutine DET (Continued) 
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from the input are used to set flags for determining possible and 
probable faulty components as described earlier. 

The program is listed showing the nearest neighbor matching 
procedure. 



D. SUBROUTINE LSQ 

Subroutine LSQ was written to determine which library signatures 
most closely matched the input signature. LSQ is similar in instruction 
flow to subroutine DET and therefore a flow chart is not presented. 

In order to determine a match, LSQ uses a method similar to the 
nearest neighbor concept. However, if one considers the signature to 
be a vector with each test frequency a dimension, LSQ does not weigh 
every dimension equally. If one dimension can be strongly affected by 
random variations in component values that dimension is not considered 
equally with the other dimensions. Thus, the distance between the input 
vector and the library vectors becomes: 

l *,<x f - V 2 



where 



t = b-jX-| + b 2 X 2 + b 3 X 3 + •• 

= input vector 

Y = C 1 Y 1 + c 2 Y 2 + c 3 Y 3 + " 
= library vector 



• th 



a. = weighting factor for the i dimension (or 
frequency) 



a. < 1 
1 — 
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LSQ operates on only one signature at a time and selects the proper 
library signature to match the input. The output is a list of the pos- 
sible signatures and the component variation defined by each one. 

One may choose to set each a . = 1 . This would be the case if the 
quantization process was sufficient to take care of the random variation 
problem. 
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VIII. SUGGESTED FURTHER STUDIES 



A. ALTERNATE APPROACH TO SIGNATURE TABLE 

Table VIII-1 lists the calculated response at each test frequency 
for the primary test circuit. Table VIII -2 lists the range of the 
response at each frequency and the maximum deviation from nominal response. 
Zero output is not included. 

Test Frequency (Hz) Response 





V V 1 


Vh 


.00048 


-2.28 Db 


12.86 Ohms 


.027 


-6.77 Db 


3.22 Ohms 


.055 


-10.85 Db 


14.32 Ohms 


.299 


r56.8 Db 


3.01 Ohms 


2.99 


-120.3 Db 


4.34 Ohms 


49.97 


-215.5 Db 


53.46 Ohms 



TABLE VIII-1. Nominal Response at Test Frequencies 



Test Frequency (Hz) 


Response Range 




Maximum Deviation 




Vg/V^db) 


V 1 


1 (fi ) 


V 2 /V 1 (db) 


V^l^n) 


.00048 


.0028 


•29.91 


3.0 


66.77 


-27.63 


+53.91 


.027 


15.58 


-46.16 


.246 


300.2 


-39.39 


+296.98 


.055 


-4.36 


-55.7 


3.0 


302. 


' -44.85 


+287.68 


.299 


-20.94 


-97.14 


.23 


300. 


-40.34 


+296.99 


2.99 


-52.28 


-160.41 


3. 


319.6 


+68.02 


+315.26 


49.97 


-98.48 


-255.5 


3. 


5338. 


+117.02 


+5284.54 



TABLE VIII -2 . Response Range and Maximum Deviation 

One can see that the responses do not cover the same range of values nor 

are the maximum deviation from the nominal response the same. From the nature 

of the network functions it can be assumed that the responses at. different 

test frequencies are not independent for a given circuit condition. Due to 

these situations it would seem logical to weigh equal changes at different 
frequencies differently. 
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The setting of quantization levels and subsequent signature 
generation is one way to accomplish this weighing. The setting of 
quantization levels based on the selected situations tends to weigh 
responses at some test frequencies less than others. However,, this 
is purely an "educated guess" method and has no analytical foundation. 

An alternate scheme was studied briefly and is suggested as the 
subject of further study. 

If one considers each test frequency as one of the dimensions in 
an n-dimensional space then each response would be a point in the n- 
dimensional space. Instead of limiting calculations to a few selected 
circuit conditions one could generate many points representing many 
circuit conditions. One could then have a large enough sample of data 
to do some statistical operations. 

In order to remove dependence the normal procedure is to multiply 
each vector, representing a point, by the inverse of the co-variance 
matrix. The result is the collection of points in a "transformed space" 
where dependency is accounted for. 

The problem is the construction of the co-variance matrix.. However, 
since the data base is now sufficiently large, an estimate of the. matrix 
can be calculated from the data. 

With a circuit under test the response would constitute a vector 
which would be multiplied by the inverse of the covariance matrix. At 
this point there are several options for fault identification. One can 
make a comparison with the stored data and make the fault identification 
using the nearest neighbor technique. An alternate approach would be to 
use an average of the nearest m neighbors when m is some positive integer. 
There are other possible procedures for matching. 
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If one is only interested in which component is faulty it might 
also be of interest to investigate any possible grouping of points which 
would be useful as in pattern recognition. The unknown could then be 
assigned to a group according to its position. Reference 5 details the 
necessary pattern recognition information. 

The above described procedure would require a computer with sufficient 
storage to process and store the information pertaining to the data base 
and covariance matrix. However such a computer is readily available for 
a nominal cost. 

B. EXTENSION OF PROCEDURE 

It soon becomes obvious that for even moderate sized circuits the 
procedure would be very difficult to use. For example, a 25-element 
circuit with 6 possible element value perturbations would yield 150 
signatures. For ease in testing the number of test frequencies should 
be as small as possible. However, with 150 signatures considerable 
signature duplication is probable unless the number of test frequencies 
is large. 

Reference 3 shows how a ladder network may be divided into subnet- 
works such that the voltage gain of a given subnetwork is independent of 
the voltage gain of all preceding networks. Therefore, testing by add- 
ing subnetworks one at a time can yield the fault component. 

This procedure uses only the voltage gain but does show one approach 
to handing large circuits. In any case, a large network would have to 
be divided into subnetworks by some scheme in order to reduce the problem 
to a manageable size. 
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IX. CONCLUSIONS 



The majority of the conclusions made as a result of this study were 
based on a comparison of the results is section VI and results published 
by other investigators. 

Most investigators reported approximately 75% isolation using one 
network function with no allowance for variation-within-tolerance of the 
-nonfaulty elements. One investigator [2] reported 78%-98% success. How- 
ever, the results were not sufficiently documented to judge the merit of 
the claim. 

The results of this study show that an additional error rate of 
17%-42% can be introduced by the within tolerance variations if the 
conventional procedures are followed. However, if sensitivity functions 
are used to aid in selecting test frequencies the error rate due to other 
element variations drops to 1 1 %-31 % for one network function. Further- 
more, using two network functions fault isolation is successful for up 
to 94% of random single fault conditions. 

The high degree of successful fault isolation clearly points out the 
advisability of using at least two network functions for testing. The 
two functions used for this study were V^/V-j and / 1^ . However, the 
designed purpose of the circuit may dictate a different choice of 
functions to be used. 

The results clearly show the improved accuracy of isolation using 
sensitivity functions as a guide in test frequency selection. In addition 
to improved isolation this method also eliminates the possibility of 
unknowingly picking test frequencies which yield the same information. 
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COMPUTER OUTPUT 

VAR V2/V1 Vl/Il 
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Figure 25. Signatures from Basic Procedure 
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Figure 26. Signatures using Sensitivity Procedure 
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Figure 27. Signatures for Two Component Variations 
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Figure 27 (Continued) 
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cn cr cr o’ CD 



COMPUTER PROGRAM LISTINGS 



REAL** WM>=GA(320) 

real MAG, N m A0,NMA0,( 2,320) 

DIMENSION 0 MC G A ( 5 ,5 ) , Ed 0 ) ,0( « ) , 0 C 5) ,° ( c ) , <= ( 5) ,UN( 5) , 
1 Y C 5) , FI ( 5) , "U ( 1 0 ) , MAG< A ,1 0 ) , MMAG(6 ,10, ?0) , r>X ( *5 ) , Y X ( 5 ) , 
2Z( 5) , X M A X ( 20 ) , XM I M ( 20) , INMAG (6, 10, 20) , 0M4G< 2,10) 

COMMON NMiG 

OATA OMEGA /? c *0. / , 0/ * *0. / , R / 5*0 . / , F I /T*0 ../ ,F/5*0./ 

OATA 01/10*0./ 

READ (5,110) PACO 
E ( 1 ) = 3 . 

F( 2 ) = . 1 7 
E( 3 ) =1 .000001 
E( A. ) =1 0.00001 
E( 5) =9 . 

F( 6>=A. 

DO 3* 1=1,6 
DO 3 P K = 1 , 10 

GO TO ( 1,2, 3, A, 5, 6, 7, 3, 9,10), K 

1 TEMP = C ( I ) 

GO TO 11 

2 F( I ) = 1 .L*T?mp 
GO to ii 

3 FA I ) = 0.9* T FMP 
GO TO 11 

4 F< I )=0. 

GO TO 11 

5 E( I ) =0 . c *T = md 
GO to 11 

6 FA I )=0.**T=”P 
GO TO 11 

7 F( I ) = 1 .2* t pmp 
GO TO 11 

* E( 1 ) = 1 .5*TFm P 

GO TO 11 

9 E< I ) = 100.*TFM£> 

GO in 1 1 

10 F( I ) =TFmp 
GO TO 3 * 

11 UN ( 1 ) = p ( A ) 

D( 1 ) = F ( 1 ) + F( A) 

D( 2) = F( 2 ) +P( 5) +E( 1 ) + c ( A)*( F( 3 ) +E(6 ) ) 

Q(3)= c (5)*(-(1 )*“(3)+E(A)*F(6) )+- c (4U*E(2) <: (E(3)+F(6) ) 

D(A)=^(5)*(F(l)^ c: (A)*c(3)*F(6)+E(2)^F(3)) 

D(5)=F(2)*E(3)*p(A)*p(5)*F(6) 

Y( 1 ) =1 . 

Y(2)= c (A)*(-(*) * c (6) ) 

Y ( 3) = ^( 3 ) * c ( c ) 

Y(4) =F ( 3 ) * E( V)* c ( 5) * c ( 6) 

IF ( K .GT. I.or. I ,r,T.i , GO TO 36 
CC WR ITF ( * , 100 ) UN( 1 ) 

CC WR I Tr (*,101) 

CC WRITE (*,102) (0(N) ,N = l ,5) 

CC WRITF (p,10?) ( 0 ( N ) , N = 1 , 5 ) 

CC. WRITP (0,10A) 

CC WRIT- (3,10*) (Y(N),N=1,A) 

C DELETE starting here for sensitivity function use 

00 1 11 1 LM = 1 , 5 
DX ( L M ) = p ( L M ) 

1111 Y X ( L M ) = Y ( L M ) 

CALL PROD (0x, c ,0,FI,P CL , I R , I E R ) 

IF ( TFR. e O.O) GO TO 12 
CC WRITE ( R , 1 0 5 ) I=* 

CC 12 WRITE (*,106) ( 0 ( N ) , E I ( N ) , N = 1 , 6 ) 

12 CONTINUF 

CALL PROP (YX,4,R,c,POM t IP , IF p) 

IF ( IER . c Q.O ) GO to 13 
CC WRIT*- (?,10 c ) I-f 

CC 13 WRITF (*,106) (* (N) ,F(M) t N = l ,3) 

13 CONTINUE 

DO 16 M = 1 , A 

IF ( ABS(EKM) ) .LF.1F-10) GO TO 14 
OMEGA ( 1 , " ) = A B S ( E I ( M ) ) 
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G o to 115 

14 OMEGA ( 1.,M )=ABS( Q( M) ) 

15 CONTINUE 

00 l 7 M= 1 , 3 

IF <£BS< F< w) ) .LF.ie-10) GO tq 16 
OMEGA (2,M)=ABS( F( M) ) 

GO TO 17 

16 OMEGA (2,M)=AB5(Q(M) ) 

17 CONTINUE 
00 IB 

00 IB N = l,4 
00 IB N0=l,4 

IF IOMEGA! w,N) .EO.OMCGt ( M,NO) .AMO.N.NP.NOI GO TO 75 
GO TO IB 

75 OMEG A ( m , n ) = 0 . 

IB CONTIMUF 

DO 21 J = 1 » 4 
00 20 L = J ,4 

IF I OMEGA! 1 , J) .G T . OMEGA! 1 ,L) ) GO TO 19 
GO TO 20 

19 OMEGA ( 1 » 5 ) = 0 M EG A { 1 , J) 

OMEGA ( 1 f ,1 ) = 0 m-g A (1,U 
OMEGA ( 1 » L ) =0 VEG A (1,5) 

20 CONTINUF 

21 CONT I N U c 

DO 22 L = 1 , 4 

IF (OMEGA! 1,1) . c Q.O. ) GO Tn 22 
l P = L + 1 

OMEGA (3, r >) = (0M^GA(l,LP)-OMt:r,A(l t L) )/2. 

OMEGA ( 3 ,l)= nM EGA ( 1 , l ) -over, * (i,E) 

OMEGA <3,2)=0MFGA ( 1 ,L )+ r ' MC GA (3,5) 

IF (OMEGA ( 3 , 1 ) . LT . 0. ) WGi (3,ll=G M rGA(l,L)/2. 

L P 2 = L + 2 

IF (LP2.GT.4) GO TO ?3 

OMEGA ( 3 , K ) = (nMcr,A ( 1 ,LP2) -riMC GA (1 ,IP») /2. 

OMFGA ( 3, 3 ) = 0 M E G A (1 ,LP)+°MEGA (?,^) 

OMEGA ( 3 , 4 ) = 0 M C GA ( 1 , L P 2 ) + 0 M EG A ( 3 , 5 ) 

IF (LP2.E0.4) GO TQ 
LP 3 = l_ + 3 

OMFGA ( 3 , 5 ) = ( OMcr.A ( 1 , l P 3 ) -OMEG A (I,t_P2))/2. 

OMEGA ( 3,4) = OM c G A ( 1 , 1 P 2 ) ♦- = 0 A (3,5) 

OMEGA (3,5) = 0 MC GA( 1 , L P 3 ) + n E G A (3,5) 

GO TO 23 

22 CONTINUF 

23 00 26 J = 1 , 3 
00 2 5 L = J , 3 

IF (OMEGA ( 2 , J ) . G T . 0 m E G A (2,L)) GO TO 24 
GO TO 25 

24 OMEGA ( 2 , 5 ) =0 M EG A ( 2 , J ) 

OMFGA (2,.))=DMTGA (2,0 
OMEGA (2,L)=0MEGA (2,5) 

25 CONTINUE 

26 C0NTINU c 

DO 2 T L=I ,3 

IF (OMFGA (2,L).EO.O.) GO TO 2 7 
l P = L + 1 

OMEGA ( 6 , 5 ) = ( o M EGA ( 2 , LP ) -ovpr.A {2,L»)/2. 

OMEGA (4,1) = OM^O A ( 2 ,L)- n Mcr 7 a (4,5) 

. OMEGA (4,2) = OM-G A (2,L)+Hm=GA ( 6 , 5 ) 

OMFGA (4,3 ) = Omeg A (2,L°)+0 V "GA (4,=;) 

IF (C M EGA(4,1) .LT.O. ) OMEGA ( 4 , 1) = n veG A ( 2 , L ) / 2 . 

L p 2 = L + 2 

IF (LP2.GT.3) GO 2« 

OMFGA (4,5)=(0 M EGA ( 2 , L D 2 ) - 0 M FG A (?,[Dn/2. 

OMEGA ( 4,3) =QM-GA (?,L D )*0”EGA (4,") 

OMEGA ( 6 , 4 ) = 0 M£ GA ( 2 » LP2 ) + OMEGA (4,5) 

GO TO 2° 

27 CONTINUE 
2B M2=0 

00 29 M= 1 , 4 
DO 29 m 1 = 3 , 4 
M2 = M 2+ 1 



29 01 (M? ) rOMFGS ( Ml , M) 

01( v 2+l)=Q M rGA(3,5) 

DO ? 91 11=1,10 
00 2°! I vf = I L , 1 0 
IF ( 01 ( I J ) .f=0.0. ) 00 TO 291 
0X = 01 ( !L ) /r>\ ( i j) 

IP ( ox ,M = , 1 . .AND. ( OX.GT. .Q.AMO.nx.LT.l .1 ) ) 
291 CONTINIJF 

00 32 *3=1,0 
DO 31 ”4=W3,a 

I c (01 ( V3 ) .LF.01 (M4I1 00 TO 31 

01(10) =P'1( M3 ) 

01 ( M3 ) =01 ( MM 
01 (MM = 01( 10) 

31 CONT ! NU C 

32 C0NTTNU c 

DO 34 M 5= 1 ,9 

IF ( 0 1 ( M 5 ) . = 0.0. > 00 TO 34 

w 6= 1 0- v p 

IF (M5.F0.1) GO TO 35 
00 33 = 

33 01 (M7) =oi ( m? + M7- 1 ) 

GO rn 3c 

34 CONTIMUF 

35 01 ( M 6 + 1 ) = 0 . 

+ 1 

C STie OF l c T I n N hhoc 

36 CALL MAG’ ( M 4 , O 1 , UN , O , Y, «4G 1 

IF (K.MS.l.OP.l.NE.l) GO TO 3 7 2 

DO 3T L m = 1 , ’ 0 

MAG ( 3 t L M ) = m a G ( 1 ,LM) 

- MAG (4,( M ) = m a G ( ? , l m ) 

37 WRITf ( 4 , 1 1 c ) m A G ( 3 , L M ) , M AG ( 4 , L m ) 



CC 


WRITr 


(8,10") 


(PI ( MP ) 


, ‘<3=1 ,M6) 




WR T TP 


( 4 , 1 1 f ) 


( 01 ( M3 ) 


, M 3 = 1 , v 6 ) 


CC 


WQ 1 T e 


( ”, 10M 






CC 


W IT- 


(8,103) 


( MAG( 1 , 


N) , M A G ( 2 , N ) , N = 1 ,M6) 


372 


HQ ^00 


I K=1 , 2 







DO 500 TL =1 , M6 

!F(MAG( I'M It ).L c .O.) GO TO 4Qo 
DM AG ( IK , IL ) = 20.*AL0010( MAG ( IK, II M 
GO TO 500 
49° d m AG( I k , I l ) =0. 

500 CONTIMUF 

DO 38 L m = 1 ,10 

IF ( M A G ( 3,LM) . L F . 0 . ) GO T n 371 
NM AG ( I , K , l m ) =m * r, ( 1 , \ M ) /mag ( 3 , L M 1 
371 I F ( M A G ( 4 , l m ) . L F . 0 . ) 90 TP 

NMAG( I , K , L m + 10 ) =MAG( 2 , LM) /MAG(4, LM) 

‘ 38 CQNTINU- 

DO 3 9 K = 1 , 1 3 

39 CALL MQM( A , 2 , 3 ,K , XMA X( ^ 1 , X M I M( K ) 1 

WRITF (4,115) ( XMAX(K7) , XMIN(K7) ,KZ=1,19) 

DO 40 K = 1 , 9 

.40 CALL QUANTA ( K , XM I m ( K ) , X MAX ( K ) , F AC 0 ) 

DO 4 1 K = 1 1 , 1 9 

41 CALI QUANTA < K , X M ] N( K) , XKAX ( K ) , F ACO) 

DO 411 T = 1 , 6 

DO 4H J=1 , 10 
DO 4 1 1 K = 1 , 3 0 

411 INMAGf I , J,K)=MMAG( I , J,K) 

DO 4 3 1 = 1 , A 
DO 4 ? .1 = 4, c 

CC 4? WRITr (M12I ( I N M AG ( I , J , K) , K=1 , * 6 ) 

42 CONT I NU C 

CC WRITf ( p , 1 1 3 ) 

43 CONTINUE 

DO 45 1=1,6 
DO 44 J = 4 , 9 
M6P=M6+1 0 

CC 44 WRIT? (3,112) ( IN M AG(I,J»K) , K=1 1 ,M 6 P) 

44 CONT I Nil? 



01 ( IL) =0. 
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cc 



WRITc (0,113) 

45 CO NT IN IJP 

WR I T c (4,112) ( ( ( INMAGU , J, K) ,K = 1, M6) , J=4 r Q) r r=l,6) 
STOP I 

100 FORMAT ( 1 5X » 1 PF1 3 . 5 ) 

101 FORMAT <6X.'V2/V1 = 

1 ' ) 

102 FORMAT ( 15X,1PE1B.5, IX, • + ' ,1PE13.5 ,1X, »S**1 +'IPFL3„5, 



11X, ' 



1 , 1 p F 13 
103 format 

1 JSCXT p I 

1 04 * FO P v AT 
1 

105 FORMAT 

106 FORMAT 
1RT PFAL 

107 FORMAT 
1 QP CQQM at 
1 0 a c ORM AT 

110 FORMAT 

111 FORMAT 

112 format 

113 FOR M AT 

115 format 

116 Format 
STOP 
END 



+ • 

1 X, • S**3 



( 1 5X , 1 P r ! 3 



1 P p 1 3 
( 6X, • 



+ 

5, 



r , f • C **-3 » 

z i » i = -- 



, 1 P F1 3 . p 
) 



» 1 X, 

i 



• ) 

S**l ♦ 



‘ IP fr T3.5,.»- s 



) 



I'M AG DA 
* flP 10 .. 3 ) 
USED ARFt*r/l°10 c l?.4) 
7.1,1 MAGNITUDE'/) 



(110) 

( ////10X , 'R p SL PART ! mag PART RFAL PART 
PART T MA G PART REAL PART I^AG PA d ti/qx 
(////IX, 'THE FR FOIJF MC 1 F S 
( //IX, 'VP /VI MAGNITUDE 
( 1.PE1 5.5, 1PF?1. 5) 

< pp .0, I c ) 

( 1X,1PD20.5 ,IP2F20. 5) 

( 1 X , 7 T 4 ) 

( 1H ) 

( 2FIA..7) 

< 5E1 5.7 ) 



SURP OUT I N c MAGI (N, OMEGA, XNUM,DFN,YMUM t MAG) 

RFAL MAG 

DIME MS 10 V OMEGA ( 10) , XNU M (5 ) ,DEN(5) , YNUMf 5 ) »MAG(4, 10) 
D n 1 I = 1 , N 
W = 0 M c G A ( I ) 

YR = YNl'M ( 1 ) - YNIJM ( 3 ) *W*W 

Y I =Y NUM ( ? ) *W-YNIJM( 4) *W**3 

DP = D F N ( 1 ) —DEN ( 3) *W*W + OFN( 5) *W**4 

DI = PEN ( 2 ) *W-DFN ( 4 ) *w**3 

YM = SOP T( YP*YP+YI "YT ) 

DM = SOR T ( DR *DR + 0 1 *0 T ) 

MAG ( 1 , I ) =XU'IM( i j /OM 
1 MAG( 2 , I ) =DM/YM 
RETURN 
END 



SUBROUTINE QUANTA ( K , QL, OH , p AC ) 
c K = TFRM IM ARP AY 
C 0L = LOWER L I M T T POR TQL . 

C 0 H= UP PER LIMIT FOR TQL. 

C FAC=FACTOR for quantization 
RF Al NMAG( 6 , 10,20) , 7 ( 5) 

COMMON NMAG 
DO 5 1=1,6 
DO 5 J = 1 , I 0 
Q(J = NM A G ( I , J,K) 

7 ( 1 ) = P AC 
7( 2 ) =F AC*3. 

7 ( 3 ) = F A C * 6 . 

Z(4)=FAC*12. 

IF (OU.LT.QL) GO TO 3 
IF (OU.GT.QH) GO TO 1 
NMAG ( I , J ,K)=0. 

GO TO 5 

1 DO 2 L = 1 » 4 

QFAC = 7 ( L ) M OH-1 . ) + QH 
I p (OU.GT.QFAC) GO TO 2 
NMAG ( I, J,K)=L 
GO TO 6 

2 NMAG ( T , J,K) =5. 

GO TO 5 
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3 DO 4 M = 1,A 
OFAC=OL-? ( M ) *( l.-OU 

IF ( QU.LT. Q C AC ) 00 TO 4 
N“AG ( I * J * K ) = M + 4 
GO TO 5 

4 NMAG ( ! t J,K)=10. 

5 CONTINUE 
RETURN 
END 



SURROUTINF NOM ( I ,J1 ,J2*K,X V 'AX*X M IN) 

RE Al N«AG( 6, 10, 20) 

COMMON NMAG 
X w AX = 1 . 

XMIN = 1 . 

DO 1 1=1,1 
DO 1 V = J 1 , J 2 

IP ( NR AG ( L , w , K ) . GT . X M A X ) X M A X= N R A G ( L * M r K 1 
IF ( MM AG( L ,M,K ) . i_ T . XVI N) X Ml N= N«AG ( L ,M ,K ) 
1 CONTINUE 
RETURN 
END 
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DIMENSION I S I 3 { l2/9/6)/JSI3(6) 

D9 1 1=1/6 
DO 1 K= 1 / 9 

1 RE AD ( 5/ 100) ( ISIG( I/K/L)/U=l/6>/( iSlG-f r+6/;K/.'L)/L=l/6) 

DO A 1=1/12 

DO A J= 1 / 9 
DO 3 K = 1 / 6 

3 JSIS(<) = ISI3( I/J/O 

WRITE (6/ 101 ) I> J/ (JBIG (K1 )/Kl= 1/6) 

CALL LSD ( I/J/JSI3) 

A CONTINJE 
ICO FORMAT ( 1212) 

101 F0R U AT(//5X/ 'THE INPUT ISJ ELEMENT *'/■• £3/.- *' VAR •* NO » » / J 2 
END 

SUBROUTINE LSD ( I 1 / J 1 / JS I G ) 

REAL I D ( 36 ) / I S 

DIMENSION lSIG(36/6)/JSIG(6)/KBlG(36//6) /.PACK 6.-) 

IF ( I 1-1 ) 1,6/ 1 

6 IF ( Jl-1 ) 1/9/1 
9 F AC ( 1 ) = 1 • 

FAC(2)=.3333 

F AC ( 2 ) ■ .3333 JFACI3) = #565 JFACIM «*A62 *FAC(.5L)*t<694 

FAC(6)=.7A5 

DO 2 J= 1 / 36 

2 READ ( 5/ 100) ( I S I 3 ( J/ I )/ 1 = 1/6) / (KSIG( J/ I! ) //I-l‘i-6 ) 

1 I F ( 1 1 -7 ) 3/7/7 

7 DO 5 U=l/36 
I D ( J ) =0 

DO A 1=1/6 

ID1*KSIG( J/ I ) -JSI3< I ) 

A I D ( J ) =ID( J)+FAC( I )*ID1**2 
5 CONTINUE 
GO TO 17 

3 DO 15 J = 1 / 36 
I D ( J ) =0 

DO 10 1=1/6 

ID 1 = I S I G ( J/ I ) - JS 1 G ( I ) 

10 ID(J) = ID(J)+FAC( I ) M 01 **2 
15 CONTINUE 
17 I S= I D ( 1 ) 

DO 20 1=2/36 
IF(ID(J)-IS) 21/21/25 

21 I F ( I 1 - 7 ) 22/8/8 

8 WR I TE ( 6/ 1 01 ) J/ ( K 3 I G ( J / lJ)/IJ = l/ 65 
GO TO 25 

22 WRITE(6/101) J/ ( I S I G ( J/ I J ) / I J= 1 / 6 ) 

25 CONTINUF 

100 FORMAT (612/ 16/512) 

101 FORMAT ( 5X/ ’ A POSSIBLE SIGNATURE IS j: ' /,T3/ 5X/,6i:2 ) 

102 FORMAT( 1H ) 

WR I TE ( 6/ 102) 

RETURN 

END 
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C MAIN PROGRAM FOR USE WITH SUBROUTINE DET (FAC NOT INCLU 
DIMENSION ISIG(12,9,6),JSIG(6),LSIG(6) 

DO 1 1=1,6 
DO 1 K= 1 , 9 

1 READ ( 2 , 100 ) ( I SIG< I ,K,L) ,1=1,6 ) , ( I SIG< 1+6, K, L) ,L=1, 6) 

DO 4 1=1,6 
DO 4 J = 1 , 9 
DO 3 K= 1,6 
JSIG(K)=ISIG( I , J,K) 

3 LSIG(K)=ISIG( 1+6, J,K) 

WRITE! 3 , 101 ) (JSIG(K) , K= 1 , 6 ) ,(LSIG(K) ,K=1,6) 

4 CALL DET ( I , J , J SI G,LSIG) 

STOP 

100 FORM AT ( 1 2 I 2 ) 

101 FORMAT!/, 5X, 'THE INPUTS ARE s ' 14, 5 12, ' - - ',612) 

END 



SUBROUTINE GET ( I , J , J S I G , LS I G ) 

DIMENSION ISIG(36,6),KSIG(36,6),JSIG(6),LSIG(6),ID(36) 
DO 3 K= 1,6 

I F ( J S I G ( K ) -5 ) 13,12,13 

12 IF(LSIG(K)-5) 13,3,13 

3 CONTINUE 

WR I T E ( 3 , 101) 

RETURN 

13 IF(I-l) 2,23,2 
23 IF(J-l) 2,33,2 
33 DO 1 K= 1,36 

1 RE AD ( 2 , 100) ( ISIG(K,L ) ,L = 1 ,6) , (KSIG(K,L) ,L = 1, 6) 

2 DO 4 K= 1,36 
I D ( K ) =0 

L D ( K ) =0 
DO 4 L= 1 , 6 

ID1=ISIG(K,L)-JSIG(L) 

ID2=KSIG(K,L)-LSIG(L) 

ID(K)=ID(K)+ID1**2 

4 LD(K )=LD( K) +ID2**2 
DO 5 L= 1 , 6 

5 N ( L ) =0 

I S= I D( 1 ) 

IT=LD( 1 ) 

DO 7 11=2,36 
I F ( L D ( II ) — I T ) 16,17,17 

16 IT=LD( II) 

17 IF! ID! II ) — IS) 6,7,7 

6 I S = I D( I 1 ) 

7 CONTINUE 

DO 10 K=1 , 36 
I X= ! K— 1 ) /6+ 1 
IF! ID(K)-IS) 35,35,36 

35 N ( I X ) = N ( I X ) + 1 

36 I F(LD!K)— IT) 37,37,10 

37 N ( I X ) = N ( IX)+1 
10 CONTINUE 

DO 1^ K= 1 , 6 
I F ( N ( K ) -2 ) 40,45,45 
40 I F ( N ( K ) - 1 ) 15,46,15 

45 WR I T E ( 3 , 102 ) K, N ( K ) 

GO TO 15 

46 WR I T E ( 3 , 10 3 ) K 
15 CONTINUE 

RETURN 

100 FORMAT (612, 16,512) 

102 FORMAT ( ' THE FAILURE IS PROBABLY COMPONENT ', I 3 , I 10) 

103 FORM AT ( ' THE FAILURE IS POSSIBLY COMPONENT «, I 3 ) 

101 FORM AT ( ' THE SIGNATURE IS NOMINAL') 

END 
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